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Abstract 

Substitutional defects in a cubic symmetry (such as a lithium defect in a KC1 host crys- 
tal) can be modeled appropriately by an eight-state-system. Usually this tunneling degree 
of freedom is approximated by a two-level-system. We investigate the observable differences 
between the two models in three contexts. First we show that the two models predict differ- 
ent relations between the temperature dependence of specific heat and static susceptibility. 
Second we demonstrate that in the presence of external forces (pressure and electric field) 
the eight-state-system shows features that cannot be understood within the framework of 
the two-level-approximation. In this context we propose an experiment for measuring the 
parameter for tunneling along the face diagonal. Finally we discuss the differences between 
the models appearing for strongly coupled pairs. Geometric selection rules and particular 
forms of asymmetry lead to clear differences between the two models. 

1 Introduction 

Quantum tunneling of substitutional defect ions in alkali halide crystals leads to particular low 
temperature properties 0. Due to their misfit in size or shape such defect ions are confined to a 
potential energy landscape with a few degenerate potential wells. At low temperatures thermally 
activated crossing over the barriers is inhibited and the defect ion passes through the barrier 
by quantum tunneling; typically at a few Kelvin hopping becomes relevant. The potential 
energy landscape in which the defect ion moves is given by the host crystal and therefore 
reflects its symmetry which for most alkali halide crystals is cubic (e.g. the fec-structure of 
potassium chloride). There are only three multi-well potentials which are consistent with this 
symmetry: twelve wells at the edges of a cube, six wells in the middle of the surfaces and 
eight wells at the corners of a cube. In all cases the edges of the cube lie along the crystal 
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axes and the multi-well structure leads to off-center-sites for the defect ion. The off-center 
position has two immediate consequences: it separates the centers of charge and leads to a 
local distortion of the crystal. Hence both an electric dipole moment and an elastic quadrupole 
moment is connected to the defect which can thus interact with lattice vibrations, external 
fields or neighboring defects. Consequently only at low defect concentrations one can describe 
the situation by isolated tunneling systems. With rising concentration pairs, triples etc. of 
defects are involved until finally one faces a complicated many body system |4|, |10| ]. 

A standard example of tunneling defects is potassium chloride doped with a small amount of 
lithium ions (KCl:Li). This system is well described by isolated defects for concentrations up 
to say 20 ppm. The minima of the system lie on the corners of a cube (d ~ 1.4 A); in the low 
temperature regime the relevant degree of freedom is thus an eight-state-system (ESS). There 
are three different matrix elements for tunneling: (i) along the edges of the cube k, (ii) along a 
face diagonal / and (hi) along a space diagonal r; edge tunneling dominates the defect spectrum 
]|] for simple geometric reasons: the edge is the shortest distance between the potential minima. 
Neglecting face and space diagonal tunneling the problem factorizes into three two-level-systems 
(TLS); this much simpler model is often used for the description of the defect. 

In this paper we want to study in how far the TLS is a good approximation for the ESS (and hence 
in how far tunneling along the face and space diagonal on one hand and the particular geometry 
of the defect on the other is negligible). We find that in most contexts the TLS approximation 
is indeed acceptable; still there are quite some experimentally observable features which cannot 
be explained by the TLS approximation. 

The plan of the paper is as follows: In section 2 a tensorial Hamiltonian for the eight-state- 
system is introduced. Physical properties such as the specific heat and the static dielectric 
susceptibility are derived and compared with the results of the two-level-approximation. In 
section 3 we discuss the coupling of the eight-state-system to static external electric or strain 
fields. We propose an experiment in order to measure directly the parameter for tunneling along 
a face diagonal. (This parameter will dominate all corrections to the TLS results.) In section 4 
strongly coupled pairs of eight-state-systems are considered. We discuss the situation by means 
of group theory and investigate in particular the relevance of different asymmetry terms for echo 
experiments. Finally section 5 gives a conclusion. 

2 Specific Heat and Static Susceptibility 
2.1 Theory 

Before going into a discussion of the ESS we would like to summarize some of the well-known 
features of a TLS. This system describes a particle in a double- well-potential depending on one 
single coordinate (instead of the three space dimensions of the ESS). The tunneling Hamiltonian 
of the TLS reads in the 'local' basis (i.e. the basis where a z represents the position 'left' or 'right' 
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of the particle) 



H = k a x . 



(1) 



The symbols c x ,a z denote the Pauli matrices. The only parameter entering the system is the 
tunneling parameter ko between the left and the right state. It is given by the WKB-formula 



k = -£ exp j- — ^2mV ^ < 0. 



(2) 



Here the energy Eq is the oscillator frequency of a well, d is the distance between the two wells, 
m is the mass of the tunneling particle and Vq is the barrier height. Since a z represents the 
position operator, all external fields couple to the product Fa z where F stands for the amplitude 
of an external field (such as static electric or strain fields and electromagnetic or acoustic waves) . 

Two characteristic features of the TLS have been studied in numerous experiments: the tem- 
perature dependence of the specific heat (Schottky-anomaly) 
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and the temperature dependence of the static susceptibility 

X s tat(T) = £-ltanh(/3/e )- 
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Let us now ask: What observable differences appear if we consider the ESS? In order to find 
an answer we derive the corresponding formulas of this system and compare them to those of 
the TLS. For the ESS three coordinates determine the wave function in the local basis { \ xyz) } 
with x,y,z = ±1; the origin of the system is located in the center of the cube and the axes point 
along its edges. The position operator now becomes a vector 
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where the symbol <8> denotes a tensor product and each factor stands for one coordinate. The 
tunneling Hamiltonian of the ESS can be decomposed into tensor products; edge tunneling along 
the x-direction for example is described by k(a x <8> 1 <8> 1). Tunneling along the face and space 
diagonal involves more than one coordinate; the corresponding Hamiltonian is given by 



H = k (1 <g> 1 <g> a x + 1 <g> a x <g> 1 + a x <g> 1 <g> 1) 

+ / (1 (g> a x (g> a x + a x (g> 1 <g> a x + a x (8> a x ® 1) 
+ r (<j x ®g x ®g x ). 



(6) 



Here k, f, r are the amplitudes for tunneling along an edge, a face diagonal and a space diago- 
nal. Let us next write down the elastic and electric momenta connected with the defect. The 
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separation of charge which arises from the off-center-position causes an electric dipole moment 

p = qr, p=— q d. (7) 
The corresponding interaction energy with an external electric field F reads 

W F = -Y,FiPi. (8) 

i 

Moreover the defect distorts the host crystal locally producing thus an elastic moment. Wiirger 
0] derived as the leading contribution a quadrupole moment 
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Remember here that rj is a component of the position operator. The corresponding interaction 
energy with a strain field is given by 

W e = -gJ2Qij £ H> ( 10 ) 

where g is a coupling constant and is the tensor of distortion produced by the strain field. 
In the static case this tensor can be derived from the exerted pressure using standard elastome- 
chanical relations ||; in the case of acoustic waves it reads = l/2(djUi + diUj), where u 
denotes the amplitude of the phonon. 

Just as for the TLS the tunneling parameters can be found with the WKB-formula (0). Assuming 
that the heights of all potential barriers are of the same order of magnitude the three tunneling 
parameters differ only because of a simple geometric reason: the distance separating the minima 
is different (the edge of a cube is smaller than a face diagonal etc.). One concludes: \k\ > \ f\ > 
|r|. Hence in a first approach one can neglect / and r. Then the Hamiltonian (||) factorizes 
and three independent TLS remain, one for each spatial direction. This is nothing else than 
the two-level-approximation keeping in mind a factor of three: one ESS (with / = r = 0) is 
equivalent to three TLS. One can then conclude: as soon as / and r are not negligible differences 
between the ESS and the TLS will arise. 

Let us go further in our analysis of the ESS. The simple structure of the problem allows the 
exact calculation of its spectrum and eigenvectors. They are given by the following scheme: 

E 3 = -3k + 3/ - r : | ) 

E 2 = -k-f + r : | + ),| -+-),| - -+) 

(11) 

Ei = k-f-r : | + +-),| + -+), | - ++) 
E = 3k + 3/ + r : I + ++) 
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with 



x,y,z=±l 



(12) 
(13) 



Here tt{ = ± denotes the parity of the wave function in the i-direction. The spectrum consists of 
four levels with a threefold degeneracy of the second and the third level. The transitions shown 
in figures la and lb are derived from the electric and acoustic interactions given above. As an 
explicit example we consider an electric field oscillating in the y-direction. Figure la gives a 
complete scheme of the possible transitions. 
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Figure la: Selection rules of the ESS for an electric field oscillating in [100]-direction together 

with the Schoenflies notation. 

Every geometric constellation has its specific selection rules; we refrain from giving a complete list 
of all these possibilities and instead propose in figure lb a schematic diagram where all possible 
transitions are shown. Note that for the electric transition one of the parities TTi changes its sign 
whereas for the acoustic transitions two parities change their signs. Hence in the accoustic case 
the product ttitt2TT3 is conserved. 
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Figure lb: Possible electric (single arrows) and acoustic (double arrows) selection rules for the 

ESS 

There are three electric transitions with almost the same frequency ~ 2k. In the acoustic case 
there are two transitions with frequency ~ 4k and two 'inside' the degenerate tripletts. (These 
transitions are visible as soon as some perturbation destroys the degeneracy, cf. section 3). The 
fact that for both couplings only one transition frequency appears is crucial for the validity of 
the two-level-approximation: by construction the TLS would not be able to reflect more than 
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one single energy. Indeed the approximation cannot account for the fact that the frequency 
induced by electric fields differs from that induced by acoustic fields. 

Let us now turn to the temperature dependence of the specific heat and the static susceptibility. 
The partition function is given by 

1 
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where rji denotes the degree of degeneracy of the eigenvalue Si and k B is the Boltzmann constant. 
The free energy of the system is then given by F = — fc^TlnZ and the specific heat reads 
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Inserting the energies given in ( |ll| ) one finds for the specific heat of the ESS: 
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(15) 
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This expression reduces to the corresponding formula of the TLS when setting f = r 

Let us now write down the susceptibility x- The response of the ESS to an external field is given 
by the commutator formula Q 

= J(b«W-P 3 (to)])8(t-io) 

with the time-dependent dipole operator p(t) = exp(i/hfft)p exp(— z /fiHt) and the Gibbs mean- 
value (•) = Tr(-exp(—(3H))/Z. In order to find the static susceptibility one has to look at the 
real part of the Laplace-transform for zero frequency. This leads to 
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(16) 
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All other elements of the susceptibility tensor are then known: the symmetry of the problem 
implies X x>x = X Vt y = X ZyZ and X x>y = X XjZ = X V)Z = 0. Again this result reduces to the 
corresponding formula for the TLS in the case of vanishing face and space diagonal parameters. 

At first sight almost no contrast between the two models emerges from the above formulas. 
Taking reasonable values for / and r (say r/f = f/k^l0 — 20%), the plots look practically the 
same. The question arises whether there is at all a simple observable feature which distinguishes 
between the TLS and the ESS. Indeed such a criterion exists. It is based on the observation that 
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for the TLS the relation k^dxX = c v holds, whereas this is not true for the ESS. One simple 
consequence is that the turning point of the susceptibility and the maximum of the specific 
heat coincide for a TLS, whereas they occur at different temperatures for the ESS. To be more 
explicit: Defining 

Tmax ■ lyrCv(T = T max ) 



- turn 
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WT 2 



cv(T 



- turn J 








(17) 



one finds for the TLS 



whereas for the ESS 



turn ; 



turn 



(18) 

as a function 



holds. It is not possible to give an explicit expression of AT = ±tum 
of / and k since transcendental equations are involved. Instead we plot in figure 2 AT as 
a function of / with parameter k. It turns out that AT depends only very weakly on k for 
0.5K < \k\/ks < 1.5K and that a linear relation AT ~ 0.85 \ f\/ks holds approximately. 



< 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
f [K] 

Fi gure 2: AT — Tf Urn — T max as a function of \f\- The different points indicate different values 
of \k\ between 0.5 and 1.5 Kelvin. The line is a linear fit with a slope of 0.85 

As the ESS is the more realistic model the deviation T m ax 7^ T urn should be visible in experi- 
ments. This will be discussed in the following. 

Note that Fiory || has shown that even for low concentrations the 'Schottky anomaly' of the 
specific heat gets distorted by the dipole-dipole-interaction between the defects, and a similiar 
effect is expected for the dielectric susceptibility. Both effects - the interaction and face tunneling 
- might interfer in a real sample. But by examing probes with different defect concentrations 
it is in principle possible to distinguish between the two effects; any change depending on a 
variation of concentration results from interactions while effects that are independent from such 
a variation arise due to face diagonal tunneling. 
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2.2 Experimental Data 



Let us now take a look at the experimental data. We will make use of the specific heat data 
from Pohl et al. ^ and the dielectric susceptibility data from Enss et al. ||. Since / will be 
the leading correction to the TLS approximation we will neglect r in the following. We then 
have two fit parameters / and k and we can look for the optimal set in order to reproduce the 
data. But indeed there is a third undetermined quantity: the defect concentration of the probe 
cannot be fixed accurately. But the concentration varies only in a very small range so we focus 
on the tunneling parameters. All parameter sets consistent with experimental data (obtained 
by fitting the relations 15 and 16 to the data) indeed confirm the relation T max < Tt urn . 

For a further comparison of the two models we will proceed as follows. The TLS is completely 
specified by the parameter ko whereas the ESS is specified by the pair (k,f). Now one can 
ask: What relation between the ESS parameter pair (k, f) and the single TLS parameter ko 
must hold in order to produce either the same maximum of c v or the same turning point of \ 
? The cy-maxima of the two models coincide if the equation k + 1.7/ ~ ko is fulfilled (which 
can be verified by differentiating and solving a transcendental equation). The turning points of 
X are the same if k + 2.7/ ~ k$ holds. This shows that an TLS-analysis yields indeed different 
tunneling parameters for the two experiments which is counterintuitive; after all the same degree 
of freedom produces both anomalies. Usually such effects were thought to be due to experimental 
uncertainty; but the above considerations show that the difference is a systematic consequence 
of tunneling along the face diagonal. Based on the difference between T max and Tt urn one can 
estimate an optimal parameter set. In order to do so one has to solve a set of transcendental 
equations, where T max and T turn are given by the data and / and k are the unknown variables. 
We have listed the result in table 1. 





6 Li 


7 Li 


k/k B 


-0.63JC 


-0.45ft: 


f/ka 


-0.08K 


-0.06A" 


Cy : (fc+1.7/)/fes 


-0.77K 


-0.55A" 


X : (k + 2.7 f)/k B 


-0.85K 


-0.61AT 



Table 1: Parameter set consistent with experimental data. In the lower part of the table 

effective values for a TLS are listed. 



At first sight the results for k look smaller than these reported in the literature. But keep in 
mind that most experiments so far were fitted by a TLS approach. Therefore we have also 
listed such 'effective' parameters emerging from an equivalent fit within the TLS model. These 
values are in line with those reported by other authors. The above parameter set is also in good 
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agreement with restrictions arising from the WKB-formula ([|). Then corresponding restrictions 
concern both the ratio 



] -d(V2-l)\ (19) 



k = 6XP { 2H 
and the ratio of the isotope effect: 

6 k \ d^Wo 



— - - exp < : -(V 6 m - V 7 m) > . (20) 

k [ 2n j 

The tunneling distance d ~ 1.4A and the mass m of the lithium are known; the barrier height 
Vq between all wells should be of the same order of magnitude. The ratios proposed above 
(f /k ~ 0.13 and 6 /c/ 7 /c~ 1.4) are then consistent with Vq ~ 200 K and Vo ~ 140 K respectively. 
This seems a rather reasonable value for the barrier height. So the above given parameter set 
is at least consistent with the given data. We would not like to go further than this statement, 
but only stress, that the ESS is not in contradiction with experimental data. Furthermore it 
reflects much better the microscopic picture we have in mind when talking about a substitutional 
defect. Still the TLS seems useful: it has almost the same temperature dependence of specific 
heat and static susceptibility and is of course much simpler. Nevertheless there appear some 
aspects which the TLS cannot reproduce at all. In the following we will discuss these aspects. In 
section 3 such properties arise from the tunneling parameters / and r; we will show that these 
features make it possible to measure the tunneling parameters directly. The properties described 
in section 4 are based on the particular geometric structure of the ESS; this structure leads to 
selection rules which cannot be understood in the framework of a one-dimensional model such 
as the TLS. 



3 Interaction with External Fields 

As mentioned above the defects exhibit both electric and elastic moments; they hence couple 
to the corresponding external fields. The main effect of static fields is a modification of the 
energy levels, whereas oscillating fields lead to characteristic transition rules. We will discuss 
the features in three steps: first we look at the defect under pressure and second under the 
influence of a static electric field. For both cases we look at the field dependence and the 
selection rules for transitions induced by acoustic and electromagnetic waves. Based on these 
results we finally propose an experiment which can determine the tunneling parameter / (face 
diagonal tunneling). Indeed in section 3.1 and 3.2 we will focus on situations where / becomes 
visible. 

3.1 The Defect under Pressure 

Let us begin by considering a defect under the influence of external pressure. As discussed 
before (see eq. @) the defect exhibits an elastic quadrupole moment Qij cx rirj(l — 5ij) which 
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interacts with the strain field e^-. The latter can be derived by elastomechanical equations 
from the pressure exerted on the system Q. The interaction energy (cf. sect. 2) reads W e = 
—gJ2ij Qij E ij' Since Qij has (by definition) no diagonal part, the diagonal part of the strain field 
is irrelevant for the interaction. This fact simplifies the discussion considerably: the diagonal 
part of the strain field depends on the pressure in a rather non-trivial way. 

The energy eigenvalue problem of a defect in an arbitrary strain field (i.e. uniaxial pressure in 
an arbitrary direction) is not analytically solvable. Yet some constellations with high symmetry 
can be solved by means of group theory. Such tractable cases appear for uniaxial pressure in 
[100]-, [110]- and [11 ^-direction. Let us write down the interaction hamiltonian for these cases. 
In order to do so one has to derive the strain field which arises from the exerted pressure |J. 
For the cases we consider these strain fields and the corresponding interaction energies are listed 
in table 2. 
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Table 2: Distortion field and interaction energy for different directions 

The off-diagonal elements e of the strain fields are proportional to the pressure; the proportion- 
ality factor depends on the crystal's elasticity modulus. The case of exerting pressure in the 
[100]-direction turns out to be trivial since the strain field then has no off-diagonal elements and 
the interaction energy thus vanishes. The eigenvalue problem of the other two cases have been 
discussed in a two-level-approximation (i.e. / = r = 0) by Gomez et al. ||. We want to avoid 
this approximation and focus on the case where the pressure is exerted in the [110]-direction 
(a general discussion can be found in [|l2|). In this situation the Hamiltonian of a defect under 
uniaxial pressure reads 

H = H Q -^{l®a z ®a z ). (21) 

We will not write down any formulas concerning the eigenvalues and -vectors (which can be 
found in reference Jl2|]); instead we will plot all interesting information and give a discussion. 
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Figure 3 shows the dependence of the energy levels on the pressure. The energies E3, E4, E5 and 
E6 grow linearly with pressure; the others vary quadratically for small pressure and then end up 
in a linear regime. All degeneracies are obviously lifted. With growing pressure the spectrum 
changes from four almost equidistant levels to four dubletts; the energy difference inside such 
a dublett is asymptotically of the order of /. The upper two dubletts are separated by a gap 
from the lower two; this gap varies (asymptotically) linearly with pressure whereas the distance 
of the dublett splitting approaches ~ 2k. 
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Figure 3: Energy levels of an ESS (in units of k) versus a = Age/d? which is proportional to the 
uniaxial pressure in [110]~direction 
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Figure 4: Energy levels of an ESS with uniaxial pressure in [110]-direction and the electric (left 
picture, single arrows) and acoustic (right picture, double arrows) transitions 

Figure 4 shows the possible electromagnetic and acoustic transitions; the pressure is chosen 
such that the coupling energy is roughly the same as the tunneling energy. In the left part 
of the scheme simple arrows indicate a dipole transition (microwaves) and in the right part 
double arrows denote quadrupole transitions (acoustic waves). The symbol F here stands for 
the electric field of the electromagnetic wave and e for the (symmetric) distortion field of the 
acoustic wave (with eij = 1 for row index i and column index j). The scheme shows that most 
transitions can be induced by an appropriate choice of the direction of the oscillating fields; in 
particular the frequency ~ / can be induced by a microwave field with F oscillating in the [110] 
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direction. 



3.2 The Defect in the Presence of a Static Electric Field 

Let us now consider the defect in a static electric field (cf. eq. (^)). The coupling energy 
between the field and the defect reads Wf = — J2i FiPi- Again the general situation of a defect 
in an electric field of arbitrary direction is not analytically solvable. Two possible restrictions 
allow such a solution. One of them is to neglect / and r; in this 'two-level-approximation' the 
Hamiltonian H = Hq — Ff factorizes into three two-dimensional problems. It is then easy to 
derive the spectrum and the transitions. This has been done by Gomez et al. || for special 
symmetries, but it is in fact possible for arbitrary directions of the electric field. The second 
possibility is to consider cases where the electric field points along the [100]-, [110]- and [111]— 
direction. The resulting problem then still has a high symmetry and group theory makes it 
possible to find a solution even for finite / and r. The interesting new feature emerging is the 
occurrence of transition frequencies proportional to / (neglecting r) in the [110]- and [111]— 
cases. These transitions are interesting since they allow a direct measurement of / (which will 
be discussed below). We will consider the [lll]-case here; the others can be found in reference 
& 

The Hamiltonian for a defect with an electric field in [11 Indirection reads 

#(lll) = fc(l<8>l®er a! + l<g><T a: <g)l + (7 a: <g>l®l) 

+ / (1 (g a x (8) a x + a x <g> 1 <g> a x + a x ® a x <g 1) 

(22) 

+ r (a x <g a x ® a x ) 

- q^F stat (l(gl(ga z + l(ga z (gl + a z (glg)l). 

The presence of the exterior field breaks the cubic symmetry (Oh) of the defect. In the case 
considered one is left with the so-called C3„-symmetry group (i.e. there is one threefold axes 
and three vertical reflection planes). Group theory then shows that the spectrum consists of 
four singletts and two dubletts since the eight-dimensional representation decomposes into four 
one-dimensional and two two-dimensional irreducible representations. Correspondingly it is 
possible to blockdiagonalize the Hamiltonian. These blocks can then be treated in (degenerate) 
perturbation theory since /, r <C k,pF. Again we refrain from writing down details and instead 
restrict ourselves to a discussion of the results. In figure 5 the eigenvalues are plotted against the 
static field F s t a t- In contrast to the level shift due to pressure there is no linear field dependence 
of any energy level. Neglecting / and r the spectrum conserves the 1:3:3:1 degeneracy 
scheme. The four levels simply spread with increasing field as 2^k 2 + (qdF/2) 2 . Considering 
finite / and r changes the situation. The most remarkable new feature is the splitting of the 
threefold degenerate states into a singlett (E\,E2) and a dublett (E^,Eq). This splitting is 
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proportional to the face diagonal tunneling parameter /; it is given by 

= 3(pF/V3)* 
6 ( P F/V3) 2 + k 2 ' 
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Figure 5: Energy levels of an ESS (in units of k) versus b = pF/y3\k\ which is essentially the 
static electric field in [111] -direction 



F II F 



F _LF„ 



8-v 



( £^ z - Ey Z ) 

(2e xy - 



E 4 

E 6 
E 2 

E 5 
E, 

E 3 





\ 










■ 
















/ 


\ 








\ 


/ ' 




s 


X 

/ 







E_, 

E„ 
E : 

E 5 
E, 

E, 















II 








1/ 


1 11 ■ 






























1 














irx 































Figure 6: Energy levels of an ESS with electric field in [111] -direction and the electric (left 
picture, single arrows) and acoustic (right picture, double arrows) transitions 

Figure 6 gives the transition scheme of the situation. In the left part of the scheme we again 
show the dipole transitions. Two cases are listed: one where the electric field of the microwaves 
oscillates in line with the static field and one where it oscillates in a perpendicular direction. 

Suppressing amplitudes oc \f /k\ 2 we find in the first case four transitions: the selection rules 
remain essentially unaltered and lead to transitions of the order of ~ 2k. If the oscillating field 
is perpendicular to the static field many transitions are possible: those 'over the gap' with a 
frequency growing with increasing static field as well as those taking place 'inside' a triplett. 
The frequency of the latter is asymptotically determined by /. The probability for the dipole 
transition grows with the amplitude of the static field oc F 



stat ' 



The acoustic selection rules are shown in the right part of figure 6. There is a whole variety of 
possible transitions; again we want to focus on the one with frequency /. The distortion tensor 
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£jj for such a transition has to be proportional to ai(2e xy — e xz —€ yz ) +a2(€ xz — € yz ) with arbitrary 
a%. One possible choice is a sound wave propagating along the (y — z)-direction being polarized 
in the x-direction. The transition amplitude of such a constellation is given by 
The increasing static field reduces the amplitude which vanishes for F S f a t 



k 2 -2(qdF stat /2) 2 



_ k 2 +2(qdF stat /2) 2 

\[2kqd. For coupling 



energies greater than the tunneling energy the transition probability then again increases. 



3.3 Measuring Face Diagonal Tunneling 

The results of the preceding sections show that there are different possibilities to measure the face 
diagonal tunneling parameter by choosing an appropriate geometric constellation for the exterior 
static and oscillating fields. Let us focus in the following on the situation of a static electric 
field pointing along the [11 Indirection (cf. section 3.2). The electric field of the microwave is 
assumed to oscillate perpendicular to the static field (for example in the [O,l,-l]-direction). The 
selection rules tell us that a transition with frequency ~ / (which should be about lOOmK) is 
possible. In order to have a transition rate of the order of unity, one has to apply a static field 
which leads to an interaction energy of the order of k. This is fulfilled if pF ~ k. Since the dipole 
moment of the defect is known to be p ~ 2.QD and k ~ IK (cf. Qj), the field strength should 
be F s tat ~ 10 6 — 10 7 V/m. To our knowledge such an experiment has never been performed. 



4 Strongly Coupled Pairs 

In section 2 we argued that a substitutional defect exhibits an electric as well as an elastic 
moment. These momenta lead to a coupling of the defects to external fields; moreover they 
produce an interaction between neighboring defects. Here again both are present, the quadrupole 
coupling of the elastic moments and the dipole-dipole-coupling of the electric moments. In 
potassium chloride doped with lithium the electric interaction dominates the elastic coupling by 
an order of magnitude. Hence only the dipole-dipole interaction 

W = ^ " 3M^) = { (e ie - 2 - 3(e 1 e fl )(e 2 e fl )) (24) 

with 

1 v 2 

J = (25) 

is relevant. Here pi denotes the dipole moment of defect i, R is the vector connecting the defects, 
the e 's are e = (2/d)f, £q and e are the dielectric constant of the vacuum and the potassium 
chloride matrix respectively. 

The structure of the host crystal confines the defects to discrete sites on the lattice. Hence the 
possible distances between the pairs will be a discrete set. Although this fact is irrelevant for 
the bulk of the pairs which have a distance of, say, ten lattice constants or more (the crystal 
then appears as a quasicontinuum) , it is most relevant for the strongly coupled pairs which we 
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will consider here. Geometric considerations show that on the fee lattice the nearest neighbors 
NN1 lie along the [1/2,1/2 ,0]-directions and the next nearest neighbors NN2 along the [1,0,0]- 
direction. Comparing the energy scales for neighboring pairs one finds that the interaction energy 
dominates the ESS tunneling parameters. This can be estimated by evaluating the dipole-dipole 
interaction ( f24| ) using the lattice constant a ~ 6.23A and the dipole moment p ~ 2.6D. For 
the nearest neighbors NN1 one is lead to a value of eJ ~ 710K (evaluated with R = a/v2). 
In order to determine the value of the coupling one has to specify the dielectric constant of the 
material. Yet on the atomic length scale it is not clear which value to take for this constant; 
at least one can confine it to the interval e(vacuum)= 1 < e < e(KCl)= 4.25. This tolerence 
leads to values 170 K < J^Ni/ks < 710 K. For the next nearest neighbors NN2 the interval is 
60 K < JNN2/ks < 250 K. This estimation, as rough as it is, shows that the coupling between 
neighboring defects exceeds by at least one order of magnitude the intrinsic energy scales (which 
are given by k « IK). The strongly coupled pairs hence constitute a composite degree of freedom 
which should be considered as a unit. 



Klein 11] and one of the authors |l0| have discussed defect pairs in the two- level-approximation; 
within this approximation one can predict the existence of a small frequency Ak 2 / J for strongly 
coupled pairs. For nearest or next nearest neighbors this energy should be between 5mK and 
50mK. Weis et al. |J investigated the system experimentally by measuring Rabi frequencies 
in spin echo experiments. He found a broad distribution of frequencies in the lOmK region; in 
addition there was a rather narrow distribution of Rabi frequencies VLr cx 4k J J Fp. The analysis 
of these experiments was first done by considering the strongly coupled pair as an effective TLS 
with a tunneling rate Ak 2 / J and an asymmetry A. The total energy of such an effective TLS 
then reads E = \J A 2 + (4k 2 / J)' 2 (which is the resonance energy of the external field) . Weis now 
argued that the measured frequency distribution is due to a distribution of the asymmetry A. 
Such a distribution arises from imperfections of the probe, whereas there is no plausible reason 
for a distribution of the dipole moment. 

Still the existence of only one frequency was somehow striking; obviously only one pair con- 
stellation contributed to the signal (otherwise additional Rabi frequency peaks would appear). 
In order to explain this puzzle Wiirger [|J] discussed a pair of ESS using perturbation theory. 
He showed that only the pair constellation NN2 has a frequency in the range between 5mK 
and 50mK. For all other constellations geometric selection rules prohibit transitions with these 
energies. He also showed that an asymmetry proportional to the position operator explains the 
distribution of frequencies. Yet he left open the question which role other forms of asymmetries 
(e.g. a 'quadrupole asymmetry' proportional to r%Tj which couples to internal strain fields) will 
play. 

In order to treat this problem one can discuss the NN1 and NN2 case applying group theory and 
tensor factorization. In this context all types of asymmetries are easily discussed considering 
them as small perturbations. We restrict ourselves to a brief overview of the proposed methods 
and present the results for the NN2 constellation (for details cf. fll2|). 
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First the basis { \xyz) } for one defect is expanded to { \x\X2y\y2Z\Z2) } for the pair. In this 
basis the Hamiltonian reads 

H = H + W (26) 

where Hq stands for the tunneling energy of both defects and W denotes the dipole interaction. 
This (64-dimensional) Hamiltonian indeed factorizes just as in the case of a single defect (we 
write ki for the edge tunneling rate of the defect i): 

H = ( h(a x tg> 1) + k 2 (l (8> <T X ) ~ ^(<7 2 <8> <T Z ) ) ® l4®4 ® l4®4 



2 

+l4(X)4 <S> l4®4 <8) ( A;i(a x (8) 1) + fc 2 (l ® o^) + ^-(c* <8> ^2 



+l4®4®(^l(T 1 ®l)+Ml®ff I ) + ^8^)) ®3-4®4 (27) 

J, 



We are left with three four-dimensional problems. After this factorization it is easy to derive 
the complete spectrum together with all selection rules. This spectrum consists of different 



multipletts which are separated by gaps of the order of ~ For the experiments considered the 
temperature is much lower than this energy and hence only the lowest level group is of interest. 
In the middle of figure 7 the spectrum of this lowest multiplett is shown together with the dipole 
selection rules. It consists of eight states with a degeneracy scheme 1:1:2:2:1:1. The 
energy eigenvalues are typically the sum of three roots of the form \J ( J/4) 2 + (k± ± /C2) 2 and 

In order to determine the transitions caused by an external field one has to consider an interaction 
with both defects: Wf = —F(p\ + P2)- The spatial dependence of F is neglected since the 
wavelength of F (~ 1cm) is much greater than the distance of the two defects (~ 10A). There 
is only one transition frequency (just as in the case of a single defect). But here transitions are 
only induced by electric fields F oscillating in line with the distance vector R. 

The structure of the spectrum together with the selection rules makes it again possible to talk 
about the pair as an effective TLS with tunneling rate \k\k2l J ■ In particular it is possible to 
use the well-known Rabi-formalism for TLS in order to interpret the echo experiments. In a 
way this is an a posteriori justification of the TLS approximation which had been proposed in 
the context of defect pairs [fTPf] . 

Considering now an asymmetry is somewhat more complicated as in the two-level-approximation. 
There the asymmetry is simply proportional to the one-dimensional position operator. The ESS 
is in a sense three-dimensional and one has to find a way to assign a different potential value to 
each well. This can be done by introducing three dipole terms, three quadrupole terms and an 
octupole term. Allowing different asymmetries at both defects one gets: 

V d = ]T Viirj + r^+Sviirj-r!) (28) 

i=x,y,z 

V q = £ ViMrj+ry^ + Sv^rjrj-rfr]) (29) 

i,j=x,y,z 

V = v xyz (r\r\r\ + r 2 x r 2 y r 2 z ) + bv xyz {rlr\r\ - r 2 x r 2 y r 2 z ), (30) 
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with the mean asymmetry v = (v\ + V2)/2 and the difference 5v = (vi —V2)/2. The existence of 
dipole and quadrupole asymmetries is plausible since it can be deduced from electric and elastic 
internal fields which couple weakly to an arbitrarily chosen defect. In contrast the octupole term 
has no such explanation and we hence neglect it in the following. The asymmetries and V q 
will be of the same order of magnitude as Ak\k2/ J ■ Treating them properly one is led to the 
following conclusions: 

(i) We find that the part of the dipole mean asymmetry which is parallel to the distance vector 
R alters the frequency to y/ (4fci k2/J) 2 + v 2 ; this result is consistent with that proposed by 
Wiirger j|. 

(ii) All other terms (such as the quadrupole terms) do not change the transition frequency; 
neither does the difference in asymmetry. 

These results are shown in figure 7. In the middle we show the lowest multiplett without 
asymmetry; on the right hand side we show the 'dipole asymmetry' changing the energy levels, 
on the left hand side we show the 'quadrupole asymmetry' simply shifting the levels without 
changing the frequencies of the allowed transitions. The experiments by Weis are hence only 
sensitive to the mean dipole asymmetry. All other forms of asymmetry do not play a significant 
role here. 

quadrupole without dipole 



asymmetry 



without 
asymmetry 



asymmetry 



Figure 7: The lowest multiplett of NN2 



5 Conclusion 

In this article we had the aim to compare the physics of a two-level-system (TLS) to that of 
an eight-state-system (ESS). Both systems are used as possible models for cubic substitutional 
defects such as the lithium atom in a KC1 host lattice. The ESS reflects much better our 
microscopic picture of the defect and it is hence interesting to know what features are left out 
of consideration when treating the defect as a TLS. 

Using a tensorial notation we reminded the reader that the ESS reduces to three TLS if tunneling 
along the face and space diagonal of the cube is neglected. Then we presented three contexts 
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where differences between the two models arise. First we looked at the temperature dependence 
of the specific heat and the static dielectric susceptibility. A simple way to summarize the 
differences of the models is to state that for the TLS the relation k 2 drX = c v holds, whereas 
this is not true for the ESS. In particular the maximum temperature of specific heat and the 
turning point of the susceptibility coincide for the TLS, whereas they are different for the ESS. 
We think that the experimental data indeed shows this tendency. The ratio f/k (face to edge 
tunneling) and the isotope effect 6 /c/ 'k estimated from the data leads to reasonable orders of 
magnitude for the barrier height (100-200K). 

Second we investigated the differences of the two models in the presence of external forces 
(uniaxial pressure and static electric field) . There appear different aspects where tunneling along 
the face and space diagonal of the cube is visible. These features make it possible to set up a 
resonance experiment where the face diagonal tunneling can be measured directly: applying the 
electric field along the [11 Indirection, there are (both acoustic and electromagnetic) transitions 
with frequencies u ~ f{pF) 2 / {(pF) 2 + k 2 ) (where p denotes the dipole moment and F the static 
electric field). Of course such features cannot be described within the framework of a TLS. 

As a third context we discussed strongly coupled defect pairs. Here again some properties arise 
that cannot be explained by a two-level-approximation. For strongly coupled pairs geometric 
selection rules inhibit transitions for some constellations (such as the nearest neighbors). In 
addition different asymmetry terms are possible for an ESS, whereas in the TLS only one such 
term appears. It is by the interplay of selection rules and the structure of the spectrum, that 
these other asymmetries are not observable in the echo experiments performed. This is why the 
situation can be described in terms of an effective TLS. 

We thus conclude that there is a number of situations where the ESS shows observable features 
which cannot be described by a TLS. Nevertheless the latter is a good approximation in all 
contexts where the geometry and the parameters for tunneling along the face and space diagonal 
are of no relevance. In these situations the TLS being the simplest model for a tunneling degree 
of freedom indeed sketches all essential features of the defect. 
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